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ABSTRACT 

An investigation ofa viscoelastic material damping effect is studied on an example ofplenum air-cushion craft 
model. A numerical investigation was conducted to determine the vertical response characteristic ofan open plenum air- 
cushion structure. The pure vertical motion ofan air-cushion structure is investigated using a non-linear mathematical 
model; this incorporates a simple model to account the hereditary deformable characteristic ofthe material. 
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1. INTRODUCTION 

The term “Viscoelastic” material has quite a broad meaning. For example, in the literature [1 to 6] there is 
a use of this term if an equation of motion includes a viscous damping term; i. e., the equation of motion is written 
in terms of current instant values of displacement, velocity, and acceleration. In this study, the equation of motion 
will include the integral term and history of strain is required for its formulation. Materials yielding such a 
constitutive relation (requiring history) are also called viscoelastic, but the term “hereditary” materials will be used 
in this study to distinguish them. 

Hereditary properties are present in any composite material [1,2 and 3]. The definition of the hereditary 
medium will give below. At this stage, we just note that the class of viscoelastic materials includes a subclass the 
hereditary materials; i. e., these two terms are not identical. 

In this study, numerical integration is applied for solution of the integro-differential equation. According to 
the correspondence (Volterra) principle, a solution of the viscoelastic problem can be obtained from a solution of 
the corresponding elastic problem be the replacement of elasticity constants by their hereditary analogs (integral 
operators). 

In this study, the use of exponential terms for relaxation kernels is utilized. 

2. MODEL OF HEREDITARY DYNAMIC SYSTEM 

Consider the disturbed motion of system near equilibrium with a view to sustaining only vertical motion. 

As to the constitutive relations, there are different models in use. Different models of viscoelastic material 
are discussed in references [4, 5, and 6]. 
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In this study, we build a model, which based on a constitutive law of the form: 

a& + <y-by + cy, (1) 
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In the references [1-6] are listed the different standard models of viscoelasticity. Therefore, the viscoelastic 
system (1) have the properties of complex viscoelastic suspension and the relation between g and y, can be written as a 
hereditary type of exponential kernel of relaxation 
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However the model of the viscoelastic system (2) with an exponential kernel (3) includes creep strain, and stress 
relaxation, but have one major weakness such as y(t ) in initial time has the final values, and do not fulfilled with 
experiment. This disadvantage easily overcome by use of weak-singularity features of the relaxation kernel (3) following 
from [3]: 
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The fact, when the model of the system is made of composite materials [4], then the relation between o and y must 
obey the law of hereditary non-linear theory viscoelasticity such: 



(>’(0 - jy 3 (o)- J R(t ~ r)\y(T) - tf(t)\lt 


(5) 


The relation (5) is fairly common because in particularity can be obtained by the standard model of a viscoelastic 
body with kernels of relaxation (3). If considered that nonlinearity is J — 0 , will be obtained the known linear relation of 
the hereditary theory. 

The equation of motion will be: 

my(t) + a = 0 (6) 

Substituting (5) into (6) is built weak singular integro-differential equation of nonlinear hereditary deformable 
system. This equation in dimensionless coordinates can be written in text form 
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where fl — -—, N = —, T(t — T) = A 0 e~^ f ~ T) ( t — T ) a ~ l . The initial conditions are 

2^t/ 0 K 

U(0) = u 09 U(0) = U 0 ( 8 ) 

The equation (7) with initial conditions (8) represents a mathematical model of the hereditary deformable system. 

3. NUMERICAL EXAMPLE 

The results of calculations of steady-state responses according to the expression of previous sections are presented 

below. 

As the example will be calculated the vertical dynamic response of plenum air cushion craft. Both theoretical and 
experimental research [9, 10, 12-14] has been performed to study the vertical motion and /or stability of various air cushion 
design configurations. The theoretical works [10] have been based on linear approaches to analyses air cushion vertical 
responses. Since air cushions are nonlinear, the application of linear analysis may be insufficient to predict fully the air 
cushion dynamic response behavior. Furthermore, any nonlinearity in air cushion response will have a direct bearing. 

In this section presents the results of a numerical investigation to determine the dynamic behavior (in vertical) of a 
simple plenum air cushion suspension system in response to steady-state disturbances. The linear and nonlinear 
viscoelastic behavior of the system is examined. 

An air-cushion craft structure is loaded by a vertical pressure load P , which is graphically presented in Figure 1 
for a mass flow rate Q ( Q in and Q out ) to cushion. 



Figure 1: A Plenum Air Cushion scheme 


From this Figure 1 geometry of the air cushion represented by h -the distance between the ground and a lower 
surface of the air cushion body and H is the average height of the air cushion. Motion parameters of the air cushion can 
be described by 


V — V 0 (1 — f~) and h = h 0 
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where y = h 0 — h (displacement of the air cushion). 
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Relation between inner volume pressure and density of air equal to pp n = COTlSt , where n is the polytrophic 
parameter of air. 

By using equation (7) with initial conditions (8) is built a mathematical model of the air cushion motion and 
stability problem. There are two unknowns: U (t) and N . Therefore, it is needed to find the values of N on which 
system is self-vibrating. 

Exact solution of the equation (7) is not possible, so the numerical simulation of this equation with initial 
conditions (8) will be done by the method, offered in [7] 

U}-TU) Bke !K [ Ujk _ ^_ k ]j (1 0) 

a k =o 

where 


n —1 
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t n — nAt, cij — A t, a 0 = a n = At /2 
j = 1, n -1 ; k = 1, j -1 
B 0 = At a /2; Bj =At a [j a -(j-l) a ]/2 
B k =At a [(k + l) a -(k-l) a ]/2 


To solve equation (10) is setup technical and kinematical characteristic of the air cushion as followed: h 0 - 
distance between the lower surface of the air cushion and ground equal [0.02-0.12] meter; average height of the air cushion 
height is 1.25 meter; n = 1.4 and pressure ratio is 1.2. Rheological and nonlinearity parameters are varied for calculation 
of elastic, viscoelastic, and linear visa versa nonlinear cases. 

Time-history of the linear dynamic response fbr N=0.02 Time histiry of the nonlmear dynamic response for N=0.02 




Figure 2: Time History of the Vertical Dynamics Responses for N=0.02 
(a) Linear Case and (b) Nonlinear Case 
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Time history of the linear respcnse for N=0,Q4 Time history of Ihe nonlinear dynamrc response for N404 




Figure 3: Time History of the Vertical Dynamics Responses for N=0.04 
(a) Linear Case and (b) Nonlinear Case 

Time history of the linear vertical response for N=3.5. Instably response Time history of the nonlinear vertical response for N=3.5 




Figure 4: Time History of the Vertical Dynamics Responses for N=3.5 
(a) Linear Case and (b) Nonlinear Case 


Tine histories of the dynamic responses for N=2.52(critrcal value) and N=3.0 



From the analysis of solutions of the equation (7) follows that jd — N > 0 in both linear and non-linear case is 
taking place damped oscillatory process. The damping speed and dissipation characteristics of the system are significantly 
dependent on the rheological parameters A 0 , /3 and a . When the system has smaller the singularity parameter Cl of the 

structure material, then damping properties of this material is higher. The self-vibrating is occurs only if jd — N < 0. 
Effect of physical non-linearity’s and rheological parameters on the critical speed A^ cr is not difficult to show, through 
computer experiment based upon an algorithm (9). 
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4. CONCLUSIONS 

Application of the numerical integration method to the hereditary deformable problem is demonstrated. In this 
study, the numerical solution in the time domain for the dynamic problem (stability and flutter problem) have been 
presented. 

The constitutive relation (stress-strain) was used in form of a hereditary law with the relation kernel represented 
by Abelian type function. 

Numerical experiments for a problem of system bending vibration under force have been conducted. Finally, an 
understanding of viscoelastic phenomena may be exploited to damping in a manner that improves stability the of the craft. 
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